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We consider a self-consistent system of Bianchi type-I (BI) gravitational field and a binary 
- - , mixture of perfect fluid and dark energy. The perfect fluid is taken to be the one obeying 

^ i the usual equation of state, i.e., p = C,e, with C, G [0, 1] whereas, the dark energy density 

O \ is considered to be either the quintessence or the Chaplygin gas. Exact solutions to the 
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corresponding Einstein equations are obtained. 
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The description of the different phases of the Universe concerning the time evolution of its 
^ . acceleration field is among the main objectives of the cosmological models. There is mounting 
evidence that the Universe at present is dominated by the so-called dark energy. Although the 
nature of the dark energy (DE) is currently unknown, it is felt that DE is non-baryonic in origin 
lUl]. It is also believed that the dark energy has large negative pressure that leads to accelerated 
bi). expansion of the Universe. 

In view of its importance in explaining the observational cosmology many authors have con- 
is^ ; sidered cosmological models with dark energy. The simplest example of dark energy is a cos- 
5_j ■ mological constant, introduced by Einstein in 1917 0]. The discovery that the expansion of the 
^ Universe is accelerating [3] has promoted the search for new types of matter that can behave like 

a cosmological constant [4, 5] by combining positive energy density and negative pressure. This 
type of matter is often called quintessence. Zlatev et al. [6] showed that "tracker field", a form of 
qiuntessence, may explain the coincidence, adding new motivation for the quintessence scenario. 

An alternative model for the dark energy density was used by Kamenshchik et al. where the 
authors suggested the use of some perfect fluid but obeying "exotic" equation of state. This type 
of matter is known as Chaplygin gas. The fate of density perturbations in a Universe dominated 
by the Chaplygin gas, which exhibit negative pressure was studied by Fabris et a l. [^] . Model with 
Chaplygin gas was also studied in the Refs. In a recent paper Kremer ill ill has modelled 

the Universe as a binary mixture whose constitutes are described by a van der Waals fluid and by a 
dark energy density. In doing so the authors considered mainly a spatially flat, homogeneous and 
isotropic Universe described by a Friedmann-Robertson-Walker (FRW) metric. 

The theoretical arguments and recent experimental data, which support the existence of an 
anisotropic phase that approaches an isotropic one, lead to consider the models of Universe with 
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anisotropic back-ground. The simplest of anisotropic models, which nevertheless rather com- 
pletely describe the anisotropic effects, are Bianchi type-I (BI) homogeneous models whose spa- 
tial sections are flat but the expansion or contraction rate is direction-dependent. Moreover, a BI 
universe falls within the general analysis of the singularity given by Belinskii et al 1 12] and evolves 
into a FRW universe 1 13] in presence of a matter obeying the equation of state p = ^ £, C < 1- 
Since the modern-day Universe is almost isotropic at large, this feature of the BI universe makes 
it a prime candidate for studying the possible effects of an anisotropy in the early Universe on 
present-day observations. In a number of papers, e.g., [14, 15], we have studied the role of a 
nonlinear spinor and/or a scalar fields in the formation of an anisotropic Universe free from initial 
singularity. It was shown that for a suitable choice of nonlinearity and the sign of A term the 
model in question allows regular solutions and the Universe becomes isotropic in the process of 
evolution. Recently Khalatnikov et al. 1 16] studied the Einstein equations for a BI Universe in the 
presence of dust, stiff matter and cosmological constant. In a recent paper JitIi the author studied 
a self-consistent system of Bianchi type-I (BI) gravitational field and a binary mixture of perfect 
fluid and dark energy given by a cosmological constant. The perfect fluid in that paper was chosen 
to be the one obeying either the usual equation of state, i.e., p = C,e, with C, G [0, 1] or a van der 
Waals equation of state. In this paper we study the evolution of an initially anisotropic Universe 
given by a BI spacetime and a bimnary mixture of a perfect fluid obeying the equation of state 
p = and a dark energy given by either a quintessence or a Chaplygin gas. 

II. BASIC EQUATIONS 

The gravitational field in our case is given by a Bianchi type I (BI) metric in the form 

ds^ = dt^-a^dx^-b^dy^-c^dz^, (2.1) 

with the metric functions a, b, c being the functions of time t only. 

The Einstein field equations for the BI space-time we write in the form 
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Here k is the Einstein gravitational constant and over-dot means differentiation with respect to t. 
The energy-momentum tensor of the source is given by 

T; = {e+p)u^u^-p5;, (2.3) 

where is the flow vector satisfying 

g;xv«^M^ = l. (2.4) 

Here e is the total energy density of a perfect fluid and/or dark energy density, while p is the 
corresponding pressure, p and £ are related by an equation of state which will be studied below in 
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detail. In a co-moving system of coordinates from (I2.3I) one finds 

= e, r/ = Ti = = -p. (2.5) 
In view of (12.51) from (12.21) one immediately obtains lfl4ll 

a{t) = Z)iTi/3exp[Xi j (2.6a) 

(2.6c) 

Here D, and X/ are some arbitrary constants obeying 

DlD2Z)3 = l, Xi+Z2+Z3=0, 

and T is a function of t defined to be 

T = aba. {1.1) 

From (12.21) for T one find 

\=^i{^-p)- (2-8) 
On the other hand the conservation law for the energy-momentum tensor gives 

£ = -^-{£+p). (2.9) 
After a little manipulations from (12.81) and ( 12.91) we find 



t = ±VCi+3KeT2, (2.10) 
with Ci being an integration constant. On the other hand rewriting (12.91 ) in the form 



(2.11) 



and taking into account that the pressure and the energy density obey a equation of state of type 
p = /(e), we conclude that e and p, hence the right hand side of the Eq. (12.81) is a function of t 
only, i.e., 

x=^-^{£-p)x = ^{x). (2.12) 

From the mechanical point of view Eq. (12.121) can be interpreted as an equation of motion of a 
single particle with unit mass under the force ^(t). Then the following first integral exists iflill : 



x=^2[S-^{x)]. (2.13) 

Here <S can be viewed as energy and ^(t) is the potential of the force ^ . Comparing the Eqs. 
(ITTOb and dTTl) one finds <^ = Ci/l and 

^ [x) = -^K£X'^ . (2.14) 
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Finally, rearranging (I2.10I) . we write the solution to the Eq. (12.81) in quadrature: 



t + to, (2.15) 



where the integration constant to can be taken to be zero, since it only gives a shift in time. 

In what follows we study the Eqs. (12.81) and (12.91) for perfect fluid and/or dark energy for 
different equations of state obeyed by the source fields. 



III. UNIVERSE AS A BINARY MIXTURE OF PERFECT FLUID AND DARK ENERGY 

In this section we thoroughly study the evolution of the BI Universe filled with perfect fluid 
and dark energy in details. Taking into account that the energy density (e) and pressure (p) in this 
case comprise those of perfect fluid and dark energy, i.e., 

e = epf + eDE, p = Ppf+PT>E 
the energy momentum tensor can be decomposed as 

= {£DE + £pf + PDE + Ppf)u^u'^ -{PDE+Ppf)8^- (3.1) 

In the above equation £de is the dark energy density, /?de its pressure. We also use the notations 
£pf and ppf to denote the energy density and the pressure of the perfect fluid, respectively. Here we 
consider the case when the perfect fluid in question obeys the following equation of state 

Ppf=Cepf- (3.2) 

Here ^ is a constant and lies in the interval C, G [0, 1]. Depending on its numerical value, C, 
describes the following types of Universes 

^ = 0, (dust Universe), (3.3a) 

C, = 1/3, (radiation Universe), (3.3b) 

C G (1/3,1), (hard Universes), (3.3c) 

^ = 1, (Zel'dovich Universe or stiff matter) . (3.3d) 

In a comoving frame the conservation law of the energy momentum tensor leads to the balance 
equation for the energy density 

t . ^ 
£de + £pf= — (cde + Cpf + PDE + Ppf) ■ (3 .4) 

The dark energy is supposed to interact with itself only and it is minimally coupled to the gravi- 
tational field. As a result the evolution equation for the energy density decouples from that of the 
perfect fluid, and from Eq. ( 13.41) we obtain two balance equations 

t . ^ 

Cde + - (cde + Pde) = 0, (3.5a) 

T 

epf+-(epf + ;?pf) = 0. (3.5b) 
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In view of the Eq. (13.21) from (I3.5bl) one easily finds 

e,f = eo/T^'+^\ Ppf = eoC/T(i+^\ (3.6) 
where £o is the integration constants. In absence of the dark energy one immediately finds 

T = a2/(i+o, (3.7) 

with C being some integration constant. As one sees from (12.61) . in absence of a A term, for 
^ < 1 the initially anisotropic Universe eventually evolves into an isotropic FRW one, whereas, 
for C, = \, i.e., in case of stiff matter the isotropization does not take place. 

In what follows we consider the case when the Universe is filled with the dark energy as well. 




FIG. 1 : View of potentials when the Universe 
is filled with perfect fluid, perfect fluid plus 
quintessence and perfect fluid plus Chaplygin 
gas, respectively. 



FIG. 2: Evolution of the BI Universe corre- 
sponding to the potentials illustrated in Fig. ^ 



In Fig. [T]we have plotted the potentials when the Universe is filled with perfect fluid, perfect 
fluid plus quintessence and perfect fluid plus Chaplygin gas, respectively. The perfect fluid is 
given by a radiation. As one sees, these types of potentials allows only infinite motion, i.e., the 
Universe expands infinitely. The Fig. El shows the evolution of the BI Universe. The introduction 
of dark energy results in accelerated expansion of the Universe. The view of acceleration has been 
illustrated in Fig. |3l 

The Fig. |4l shows the evolution of a BI Universe filled with a perfect fluid and a Chaplygin gas. 
Here "d", "r", "h" and "s" stand for dust, radiation, hard Universe and stiff matter, respectively. 
As one sees, even in case of a stiff matter the Universe expands rapid enough to evolve into an 
isotropic one. 



A. Case with a quintessence 

Let us consider the case when the dark energy is given by a quintessence. As it was mentioned 
earlier, a new type of matter, often known as quintessence, can behave like a cosmological constant 
and was constructed by combining positive energy density and negative pressure and obeys the 
equation of state 

Pq = WqCq, (3.8) 
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where the constant Wq varies between —1 and zero, i.e., Wq G [— 1, 0]. In account of (13.81) from 
(13. Sal) one finds 

eq = eoq/T(i+"^\ = Wqeoq/T(i+^^), (3.9) 

with Coq being some integration constant. 

Now the evolution equation for t (12.81) can be written as 

As it was mentioned earlier the Eq. (13.101) admits exact solution that can be written in quadrature 

as 

f dx 

J , ==? + ro. (3.11) 



/Ci+3fc(eoT(i-^)+eoqT(i-^'^)) 
Here Iq is a constant of integration that can be taken to be trivial. 

B. Case with Chaplygin gas 

Let us now consider the case when the dark energy is represented by a Chaplygin gas. We have 
already mentioned that the Chaplygin gas was suggested as an alternative model of dark energy 
with some exotic equation of state, namely 

p, = ~A/£,, (3.12) 

with A being a positive constant. In view of the Eq. (13.121) from (13. Sat one now obtains 



eoc/T2+A, pc = -A/Jeoc/T2+A, (3.13) 



with Sqc being some integration constant. 
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FIG. 5: View of energy density and corre- 
sponding pressure when the Universe is filled 
with a perfect fluid, quintessence and Chaply- 
gin gas, respectively. 




FIG. 6: View of energy density and cor- 
responding pressure when the Universe is 
a binary mixture of a perfect fluid and 
quintessence and a perfect fluid and Chaply- 
gin gas, respectively. 



Proceeding analogously as in previous case for T we now have 

T = ^ + Veoc+AT2 + A/ Veoc+AT2) . (3. 14) 

The corresponding solution in quadrature now has the forms: 

dx 



^Ci+3K:(eoT(i-^) + VeocT2+AT4) 
where the second integration constant has been taken to be zero. 



t, (3.15) 



IV. CONCLUSION 

A self-consistent system of Bl gravitational field filled with a perfect fluid and a dark energy 
has been considered. The exact solutions to the corresponding field equations are obtained. The 
inclusion of the dark energy into the system gives rise to an accelerated expansion of the model. 
As a result the initial anisotropy of the model quickly dies away. Note that the introduction of the 
dark energy does not eliminate the initial singularity. 
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